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Abstract 


Daily observations of the photospheric rotation rate using 
the Doppler effect have been made at the Stanford Solar Observatory 
since May 1976. These observations show no daily or long period 
variations in the rotation rate that exceed the observational error 
of about one percent. The average rotation rate is the same as that 
of the sunspots and the large-scale magnetic field structures. 


*Now at Universite de Liege, Institut d' As trophy sique, Cointe- 
Ougree, Belgium. 
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Measurements of the photospheric rotation rate using the 
Doppler effect have been the subject of controversy ever since they 
were first made (e.g. Plaskett 1915; DeLury 1939; Hart 1954; Howard 
and Harvey 1970; and Plaskett 1973) and have been recently reviewed 
by Howard (1978). Daily observations at the Stanford Solar 
Observatory (SSO) have been made since May 1976. Analysis of the 731 
full-disk low-resolution dopplergrams made to date yields results that 
conflict with some other recent sets of observations. The SSO 
observations show no daily or long period variations in the rotation 
rate that exceed observational error, and the average rate is about 
the same as that of the large-scale magnetic field structures and 
sunspots. These conclusions are apparently inconsistent with the well 
known results of Howard and Harvey(1970) and the more recent reports 
of Livingston and Duvall(1979) and Howard et al .(1979) . 


*Now at Universite de Liege, Institut d ' Astrophysique , Cointe-Ougree , 
Belgium. 
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&£ Seine fl ot a t ion 


The SSO instrument design is very sirailiar to the solar 
magnstograph at the Mt. Wilson 150 foot tower telescope (Howard si. 
al . 1968 ; Howard and Harvey 1970; and Howard 1976) but was designed 
for only low spatial resolution observations. Since some of the 
differences may affect rotation measurements, the relevant details 
will be described. The overall design philosophy of the SSO has been 
described by Schemer et al . ( 1977) . Preliminary rotation results and 
the modification to the instrument to allow full disk 3-arc-rainute 
magnetograms and dopplergrams are described by DuvalK 1 979 ) . 
Svalgaard et al.(19?8) made a preliminary report of the conclusions to 
be described more fully in this paper. The results presented here 
differ slightly from the two earlier reports due to recalibration of 
the position, intensity, and velocity data. The calibration methods 
used in the present work are described below. 

Position Calibration 

The SSO telescope is outlined in Figure 1 . The instrument 
consists of a vertical telescope with focal length of 6.5m feeding a 
vertical 22.8m Littrow spectrograph. The spectrograph is used in the 
5th order with a reciprocal dispersion of 12.9 mm/8 at 5250 8. The 
entrance aperture consists of a 6mm square mask, an image slicer, and 
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a 0.75mm by 100mm slit. There are several systematic effects that 
must be examined when using this entrance aperture arrangement for 
low-resolution observations. These influence the effective position 
of the aperture on the sun and the measured velocity. 

The scanning system used is similiar to most magnetograph 
instruments. A limb guider controls the orientation of the 
second-flat mirror with an analog servo system. The guider image is 
produced by an auxiliary lens with focal length of 7.0m. The sun is 
scanned by moving the guider north-south and east-west in steps of 
0.001 inch. The position angle of the solar pole relative to the 
scanning axes is computed from the date and location of the coelostat 
mirror. The magnetogram/dopplergram scan is made by stepping the 
image to fixed positions and integrating for 15-seconds at each 
location, instead of with a continuous motion. The scan is made in 
east-west lines (on the sun) stepping 90 arc-seconds between 
integrations with 1 80 arc-seconds between lines. The aperture is 180 
arc-second square and is oriented parallel to the entrance slit which 
is north-south in the room. Prior to February 15, 1978 the aperture 
was rotated to be parallel to central meridian on the image resulting 
in non-uniform illumination of some of the image slicer segments. 
Although no difference in the velocity data was found for several 
pairs of observations with the aperture oriented each way, there is 
less risk of systematic errors with the present alignment. 
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In order to measure solar rotation the conversion between 


scanner coordinates and effective heliocentric coordinates must be 
accurate. The telesoope scale is defined by the 6.97m distance 
between the guiding lens and the scanning table. The angular size of 
the solar disk is computed for the day of observation and the 
geometrical posit -on on the disk is found for each scanner position. 
If the observing aperture were only a few arc-seconds wide this would 
be sufficient to define the heliocentric position, however, with a 180 
arc-second aperture the limb darkening across the aperture will weight 
the effective position of the aperture toward the center of the disk. 
An effective position is computed for each actual aperture position 
using the average observed limb darkening function. For the aperture 
position at the limb, the correction is about a third the aperture 
size. The correction decreases rapidly away from the limb. For the 
rotation analysis described below, only the center half the area of 
the disk is used, i.e. only points within 0.75 of a radius of the 
center of the disk. The maximum position correction for the points 
used is 4 arc-seconds. If the correction were omitted from the 
rotation analysis the computed equatorial rotation rate would be 0.25 
percent too low (about 5 ra/s) . 

Intensity Calibration 

As part of normal magnetograph operation the brightness of the 
image is recorded during each integration. This is primarily used to 
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compensate the magnetic field signal for brightness fluctuations. For 
each scan the individual intensity measurements are corrected for limb 
darkening and averaged over the center half of the disk. The varianoe 
of this quantity is used as a part of a quality figure for the scan. 
The magnitude of the average intensity will depend on the zenith 
angle, the photo-multiplier tube voltage, the sky quality and the 
cleanliness of the optics. Since the phototube voltages are rarely 
changed (twice per year) this quantity can be used to monitor the 
cleanliness of the sky and optics. When it was discovered that the 
average intensity could also be used to compute a correction for the 
effect of scattered light on the velocity observations, an effort was 
made to calibrate the intensity for the phototube voltages used since 
the start of the observations. 

When the telescope is -not making a magnetogram or some special 
observation, it is used for mean solar magnetic field observations. 
There are, therefore, several measures of the intensity of integrated 
sunlight for each day. All of the intensity observations for each day 
were plotted on a log scale as a function of computed airmass. A 
typical plot is shown in Figure 2. It can be seen that the expected 
linear relationship is found. The intercept corresponds to the 
intensity that would have been observed without the intervening 
atmosphere. From a plot of the intercept value as a function of date, 
the intensity changes resulting from phototube voltage changes were 
found. The calibrated disk average Intensity from the magnetograms is 


5 


plotted in Figure 3. Most of the variation with time is due to the 
varying amount of dirt on the optics. The arrows show the times that 
the coelo3tat mirror and lenses were cleaned. 

Velocity Calibration and Corrections 

The velocity observations are made by measuring the position of 
the spectral line. The spectrograph exit aperture consists of a pair 
of prisms mounted such that each prism gets light from one wing of the 
line and directs it to a photomultiplier tube. The prism assembly is 
mounted on accurate ways and positioned by a precise measuring screw. 
A shaft encoder is used to record the position of tne prisms while a 
servo system keeps them centered on the line. One count of the 
encoder corresponds to 1/4096 of a turn of the screw with 10 turns per 
inch. Thus, for the dispersion at 5250 8, the velocity calibration is 
2.752 m/s per encoder count. This corresponds to about 0.6 micron of 
physical displacement of the prisms and is about the practical limit 
for a mechanical measuring engine. The RMS variation due to 
spectrograph seeing, electronic noise, and servo errors is typically 
lOm/s for a 15-second observation in integrated light. The velocity 
calibration is measured each day as part of a series of automatic 
checks, however, since the uncertainty in the daily measurement of 
dispersion is larger than the variations, a constant value is used for 
the data reported here. 
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Before the observed line position data ean be interpreted as 
motions of the sun, the relative motion of the observatory and the sun 
must be removed and speotrograph drifts must be romoved. The Doppler 
shifts resulting from the earth's rotation and orbital motion are 
removed for each point observed using the method described by Howard 
and HarveyO970) and Howard et al .(iQ7Q). 

Other sources of drift in the velocity data include a slow 
motion of the diffraction grating during the first hour after it is 
repositioned, changes in the index of refraction of the air in the 
spectrograph, and drifts in the relative gain of the phototubes. To 
reduce the effect of the grating drift, the start of each scan is 
delayed for about one half hour after the grating is positioned. In 
order to remove the remaining drifts from the data a pair of 
five-minute integrations are made at the center of the disk, one 
before and one after the scan. It is assumed that the drifts are 
linear in time and these two disk-center averages are used to define 
the zero line. Figure i| shows the total amount of this drift for each 
of the observations used for this study. 

Once the relative motions and drifts are removed the rotation 
rate of the photosphere can be computed. The first question one must 
address when measuring solar rotation is how to describe that 
rotation. The traditional method has been to assume that the 
anti- symmetrical part of the dopplergram can be described by a 
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relation of the form; 


P * 

0 = a + b sin fa 8 + o ain't* * other terma 

whore w is the angular velocity, B is the heliographic latitude, and 

a,b, and o are to be found from a least-squares or other fit to the 

data (e.g. Howard and Harvey 1970). We use the same general method 

with a few differences. Since tne observed quantity is a 

line-of-sight velocity and some of the systematic errors are shifts in 

apparent wavelength we choose to expross the determined coefficients 

and thus the rotation rate in m/s rather that radians/s. This 

necessitates conversion of the final result to angular velocity for 

comparisons but is more convenient for studying the various systematic 

effects. The 1 v* b red shift is assumed to be symmetrical about 

central meridian so that it will not affect the rotation rate terms 

and is not explicitly removed. The magnitude of the limb shift is 

expected to be about 30m/ s at the boundary of the 0.75R disk used for 

o 

the rotation rotation studies and will be discussed in a future paper. 
Since an error in the computed position angle of the solar pole will 
produce an apparent north-south velocity at the expense of rot cion, 
the linear anti-symmetrical part of the velocity in the north-south 
direction is also found. 

The form for differential rotation used is the same as used by 
Howard and Harvey(i970) but we make the assumption that the terms b 
and c in (1) nave the same value. There is a computational problem 
resulting from the nonorthogcmali ty of the terms used to describe 
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rotation (Duvall and Svalgaard 1978 and Sfeenflo iy?7), Tho effect 
of this it -.hat noise in the data will shift the amplitude betwoon the 
b and e terms. Since the time averaged values of b and e are about 
the same when (i) i3 used to determine rotation, it is assumed here 
that b=c . Any errors resulting from this assumption will appear as 
stable features in tho residuals and w ’ : .e examined in a later 
paper. In any case, the assumed equality of b and e has a negligible 
effect on the a coefficient. 

Since the observed quantity is the line-of-sight component of 
velocity, and since rotation results in a purely azimuthal motion, the 
correction factor cosB-oosB ‘sinL is used where B is the heliographie 
latitude of the disk center and L is heliographie longitude measured 
from disk center. Thus, the relation used to express rotation has the 
form : 

2 2 

V = (e + b«sin B*(l+sin B)), cosB#cosB o .sinL + p.sinB 
The coefficients e,b, and p in (2) are found for each 3can by the 
method of least squares with each of the data values weighted by the 
intensity corrected for limb darkening. This weighting helps insure 
that an occasional cloud will not adversely affect the results. Once 
the coefficients above have been found, the error in the position angle 
can be found from: 

p err = tan ~ 1( P /e) 

and the equatorial rotation velocity is found from: 

n r 

a = \J e + p 
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in practice p Gj , r has boon found to be email with a typical value of 
O.i degree, la addition to the terms o, b, and p found above, a 

V* » 

term a Q is found from a second fit to the data with the term b 

replaced with a constant -300m/o. The coeff iccntu a. , a, and b have 

o 

been found for all reaaonaole completed scans to date, Of the total 
975 nouns 'ihh were eliminated because Lhev were Incomplete or were 
made during ouch bad weather that the data was obviously not useful. 

In addition to the above corrections, there are two types of 
systematic effect that should be considered for the SSQ instrument. 
These are the alignment of the entrance slit and the exit slit and the 
presence of significant amounts of scattered light for the earlier 
observations. (We now keep the optics scrupulously clean so that 
scattered light is of the order of one part in a thousand). 

The combined effect of the use of an image slioer, a 
3-arc-minute aperture, and the north-south direction of the entrance 
slit results in a signal in the same sense as solar rotation if the 
exit slit is not exactly parallel to the entrance slit, Prior to 15 
February 1978 there was a small misalignment that caused such an error 
(as described by Duvall 1979). In addition to altering the observed 
rotation rate, the misalignment of the slits results in an apparent 
change in rotation with distance from disk center. The net effect can 
be clearly seen by comparing uneorrccted dopplergrams for 
observations before and after the slit readjustment. Figure 5 shows 
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two such dopplergrams. The closed contours at the east and west limb 3 
In Figure 5a are the result of the misalignment. An analysis of the 
correlation between the equatorial rotation rate and intenstiy (as 
described below) shows that the data before the slits were properly 
adjusted should be increased by 2,1 percent. This correction has been 
included in the results presented here. 

£££&afr. ja£ S sa t t sc sd Li.gh.fr. 


The Problem 

Figure 6 3hows the observed equatorial rate a Q corrected for 
all the effects described above. The general character of the data is 
a sequence of slow declines followed by sudden increases. The arrows 
show the times that the coelostat mirror was removed and cleaned or 
resurfaced. There is a clear connection between the times of cleaning 
and the times of sharp increase in observed rotation rate. The 
correlation between the observed brightness shown in Figure 3 and 
Figure 6 is also obvious. After each cleaning of the optics, the 
intensity and a Q both jump. As dust and smog accumulate on the 
mirrors and lens, both the intensity of the image and the observed 
rotation rate decrease. This correlation is to be expected since as 
the dust scatters light out of the main beam, it mixes light from the 
entire disk with light from a particular region. 
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The details of the expected effect depenc on the density of 
scattering sources, the size and height distribution of the particles, 
and the angular aperture of the instrument. In any case it is 
expected that the effect on **ne rotation rate will be proportional to 
the scattered brightness just off the limb. Because of the limited 
angular aperture the effective intensity of the scattered light will 
be proportional to the density of scatterers while the direct 
intensity will decrease exponentially with the number of scattOiers. 
Therefore, in the analysis below, the natural logarithm of the 
intensity is used for comparison to the rotation rate. The effect of 
scattered light from the instrument and the sky on rotation 
measurements has been known for a long time (eg. Plaskett, 1915 ). 

The Correction 

Since there is an effect on the observed velocities due to 
scattered light, observations should only be made with clear skies and 
clean optics. In the past at the SSO this condition was not always 
met but since the observations made in less than optimum conditions 
contain a significant amount of information about solar motions, we 
have attempted to remove as much of the effect of scattered light as 
possible. This is done by finding an empirical relation between the 
intensity I and a Q . We expect: 

T T - T 

1 = 1 e 
o 

where T is the optical thickness. The scattered light intensity S off 
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the limb (as a fraction of the di3k center intensl.ty) should be 
proportional to T, so we can write: 

* - 1 0 e " k,s 00 
The decrease in observed rotational velocity should be proportional to 


S so we can write: 


V obs “ V actual -< * S 


then combining (3) and (4) we expect: 

V . = V . , -oyk (Ln I - Ln I) 

obs actual o v 3 

Since V , and Ln I are observed. 'V/k and Ln I are needed to find 
obs o 

V actual* figure 7 shows a scatterplot of V Qbg and Ln 1. The line is 
the least-squares linear regression line V = 134 Ln I + 885. A 
separate analysis for data before and after February 15, 1978 yielded 
the 2.1 percent correction to compensate for the slit adjustment 
described above. 


To complete the correction we need to know I , that is the 
intensity (corrected for phototube changes) that would have been found 
with clean sky and optics. Since August 1978 the scattered light 
intensity just off the limb has been measured at the start of each 
scan. Figure 8a shows Ln I vs S. The intercept of the best fit line 
yields Ln I = 8.56, Thus we have: 

V corrected = ^obs + < 8 -56 - I) 

Figure 8b shows the correlation between the observed rotation velocity 
and the measured scattered light. The roughly linear relationship 
here justifies the form of eqn. (4). The slope of the line in figure 




«b is the term oc in eqn (4) and similarly Figure 8 a provides directly 
a value for k in eqn (3). These combine to yield *</k=1 19 which is in 
close agreement to the value for <*/k derived from the direct 
comparison of intensity and rotation in Figure 7- Due to 
uncertainties in asser.-'.-nv the relative errors in the computed 
rotation rate and the ob* «-r*. ed scattered light, the slopes of the 
lines in Figure 8 a and 8 b have poorly determined errors. The expected 
maximum value for the slope in Figure 8 b is about 20 ra/s per percent 
scattered light. Inspection of Figure 8 shows that a smaller 
value of the slopes of the lines is consistent with the data, with the 
expected value of 20 m/s per percent scattered light, and with the 
ratio «/k found in Figure 7 , 

Discussion 


The Variability of Solar Rotation 

Figure 9c shows the equatorial rotation rate a Q as corrected 
with relation ( 6 ). Figure 9a is the logarithm of the brightness data 
from Figure 3 and Figure 9b is the data from Figure 6 shown here for 
comparison. Several statements concerning the nature of solar 
rotation can be made by examining this data. The first and most 
certain is that the day-to-day variation in the observed value is 
small with an rms variation from the mean of 26 m/s or about 1.3 
percent. During intervals of good weather, as indicated by small 
variance in the intensity, it can be seen that the solar contribution 
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to the observed velocity variations i3 probably smaller that 20 ra/s. 
These relatively small upper limits for the short term fluctuations in 
solar rotation are noticably smaller than the values previously 
reported for doppler shift type observations. 

Examination of Figure 9 also shows that there are no 
identifiable seasonal or yearly variations during the 3 years of 
available data that exceed 20 m/s. This can be confirmed from Figure 
10 which is the result of a 108-day low-pass filtering of Figure 9c. 
Considering the uncertainties in the corrections applied, there is no 
justification for concluding that even these small variations are of 
solar origin. There may, of course, be long term variations in the 
rotation rate that are masked by uncertainties in the corrections, but 
variations on the order of 1 percent or larger seem to be excluded by 
this analysis. 


Actual Value of Solar Rotation 

Newton and Nunn (1951) reported the average rotation rate of 
long-lived recurrent sunspots. Using their form for differential 
rotation, the equatorial rate from sunspot observations is 2020 m/s. 
An inspection of Figure 9 shows that, to the certainty of the 
corrections, the photospheric material also rotates at this rate. 
Considering the possible uncertainty in the corrections it is 
worthwhile to examine seperately the data since May 18, 1979. While 
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this only represents 60 observations, the scattered light has been low 

enough that no correction is necessary. Figure 11 shows this data. 

There has been no change in the procedure during this time except to 

keep the optics clean to exclude an effect from scattered light. The 

average value of the uncorrected rate a for this interval is 2016 + 

o ” 

13tn/s. In order to be fully confident of the constancy of this value, 
an additional year or more of observations must be made, but it is 
safe to 3ay that in June 1979 the photospheric material rotated at the 
Newton and Nunn rate. 

Conclusions 

Several points need to be reviewed. At the SSO, measurements 
of the solar rotation rate have been made since May 1976. In the 
earlier observations, the measured velocities were contaminated by 
scattered light. This is exibited by an obvious correlation between 
the observed intensity and velocity data. When the effect of 
scattered light is removed by empirical methods it is found that: 

A) The day-to-day variations in the equatorial 
rotation are at most 20 m/s rms about the mean, 

B) There are no long term trends exceeding 1 percent 
that can be associated with the sun rather than 
the instrument, 

C) The average rate for the photosphere is most 
likely about the same as the Newton and Nunn rate 
for sunspots. 
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It was also found that the fewer the number of terms used to express 
solar rotation the smaller the variance in the rate determined. When 
the term for differential rotation is also fit to the data the 
rotation is found to be: 

V s 2030 - 278 sin 2 B - 278 sin^ B m/s 
(at one solar radius) which is the same as: 

<* = 2.917 - 0.40 sin 2 B - 0.40 sin 4 B radians/s. 

These conclusions are different from the commonly held belief 
that the photosphere in general rotates slower than large scale 
magnetic field structures and that the rate shows substantial 
variations. It is unlikely that magnetograph-type observations made 
at other observatories suffer from the specific errors that formerly 
existed at the SSO (e.g. scattered light) but there may be some as yet 
unidentified sources of instrumental noise. If there is some noise 
which is not purely random then the observations may not be in 
conflict, but only the interpretations. 
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Figure 1 , A schematic diagram of the Stanford solar 
telescope. Note the separate light path for the 
guider image and the observing image. The actual 
distance between the guider lens and the guidor 
photodiodos is used to define the telescope 
scale. The spectrograph head is at ground level 
with the Littrow lens and diffraction grating at 
the bottom of a 23m pit. 


Figure 2. The change in observed image brightness as a 
function of airmass on a typical day. The 
intercept for zero airmass provides the 
calibration of the instrument intensity signal 
for that day. 


Figure 3. Observed image brightness for each magnetogram 
through July, 1979. The arrows show the dates on 
which the optics were cleaned. Intensity changes 
due to other instrument changes have been 
removed. 


Figure 4. Drift in the velocity signal ascribed to 
spectrograph and electronics drifts. The total 
amount of the drift between the reference 
integrations is shown. The average magnitude of 
the drift was 43 m/'s per hour. The algorithm for 
positioning the grating was changed in October 
1976 to reduce the residual grating motion. 


Figure 5. Solar velocity contours as observed. Dashed 
contours are motion toward the observer, solid 
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away. The contour interval i3 200 m/s, A solid 
body rotation would result in straight vortical 
eontoura, Data for two days is shown: a) (15 

February 1978) observed with an entrance slit 
alignment error that affects velocity 
measurements near the limb, b) 06 March 1978) 
typical for observations with the correct 
alignment . 

Figure 6. Solar equatorial rotation rate a for each 

0 

magnetogram scan included in this analysis. 
There are tick marks at the first of each month. 
The arrows show the dates of cleaning and 
realuminizing the mirrors. 

Figure 7 . A scatterpiot of the equatorial rotation rate a 

0 

as a function of the logarithm of the observed 

image brightness ln(I). The best fit line shown 

is used to determine the error in a which is due 

0 

to the correlation with intensity. 

Figure 8. Scatterplots of a) ln(I) and b) a with the 

0 

observed scattered light S for those observations 
for which S was measured. S is the brightness in 
the 3-arc-minute aperture centered about 2 
arc-minutes off the limb relative to the disk 
center brightness. The intercept in 8a shows the 
brightness with no scattered light I , while the 
slope of the best fit lines in 8a and 8b verify 
the correction factor derived in Figure 7. 
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Figure 9. a) The observed image brightness (from Figure 3), 
b) the observed rotation rate a o (from Figure 6), 
and o) the corrected equatorial rotation rato for 
all included scans. The correction wa3 applied 
separately for each 3can thus somo of the 
variance in brightness which is not correlated 
with scattered light tends to increase the 
variance in the corrected rotation rato. The 
estimate for day-to-day variations in a Q is made 
from Figure 9b disregarding the drift resulting 
from the correlation with 9a. 

Figure 10. Slow drifts in the oorrec'.ed equatorial rate. A 
108-day low-pass filter has been applied. The 
trends in this figure are probably due to the 
Incomplete removal of instrumental problems. 


Figure 11. The observed (not-corrected) equatorial rate for 
the most recent data (clean optics). The 
horizontal line is drawn at the Newton and Nunn 
rate of 2020 m/s. The simple mean of the data 
shown is 2016 m/s with a standard deviation of 13 
m/s. 
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